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Abstract 

The low-lying even-parity states of the nucleon are explored in lattice QCD using the PACS-CS collaboration 2+1- 
flavor dynamical-QCD gauge-field configurations made available through the International Lattice Datagrid (ILDG). 
The established correlation-matrix approach is used, in which various fermion source and sink smearings are utilized 
to provide an effective basis of interpolating fields to span the space of low-lying energy eigenstates. Of particular 
interest is the nature of the first excited state of the nucleon, the N\ Roper resonance of P\ \ pion-nucleon scattering. 
The Roper state of the present analysis approaches the physical mass, displaying significant chiral curvature at the 
lightest quark mass. These full QCD results, providing the world's first insight into the nucleon mass spectrum in the 
light-quark regime, are significantly different from those of quenched QCD and provide interesting insights into the 
dynamics of QCD. 
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1. Introduction 

The first positive parity resonance of the nucleon, the 
(1440) or Roper resonance, has been the subject of 
extensive interest since its discovery in 1964 [1]. This 
P-wave isospin-1/2 spin-1/2 (Pn) pion-nucleon reso- 
nance has held the curiosity and imagination of the nu- 
clear and particle physics community due to its surpris- 
ingly low mass. For example, in constituent quark mod- 
els the lowest-lying odd-parity state occurs below the 
Pn state whereas in Nature, the negative par- 

ity A^2 _ (1535) Si i state is almost 100 MeV above the 
Roper resonance. 

This phenomenon has led to wide speculation on the 
possible exotic nature of the Roper resonance. For ex- 
ample, the Roper resonance has been described as a hy- 
brid baryon state with explicitly excited gluon field con- 
figurations yj, |5[], as a breathing mode of the ground 
state |0] or a state which can be described in terms of a 
five quark (meson-baryon) state 0]. 

The elusive nature of this low-lying resonance is not 
constrained to model calculations alone. There have 
been several investigations of the low-lying nucleon 
spectrum using the first-principles approach of lattice 
field theory. 



The lattice approach to Quantum Chromodynamics 
(QCD) provides a non-perturbative tool to explore the 
properties of hadrons from the first principles of this 
fundamental quantum field theory. Numerical simula- 
tions of QCD on a space-time lattice with the light up, 
down and strange dynamical-quark masses similar to 
those of Nature are now possible [8]. As such, some 
long-standing problems in nuclear-particle physics are 
now being resolved. For example, the ground-state 
hadron spectrum is now well understood J9J] . 

However, gaining knowledge of the excited- 
state spectrum presents additional challenges. The 
Euclidean-time correlation function provides access 
to a tower of energy eigenstates in the form of a sum 
of decaying exponentials with the masses of the states 
in the exponents. The ground state mass, being the 
lowest energy state, has the slowest decay rate, and 
is obtained through the analysis of the large-time 
behaviour of this function. However, the excited states 
appear in the sub-leading exponentials. Extracting 
excited state masses from these exponents is intricate 
as the correlation functions decay quickly and the 
signal to noise ratio deteriorates rapidly. In addition, 
the spectrum is composed of both single-particle states 
and multiple-particle states interacting and mixing in 
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the finite physical volume of the lattice. Understanding 
the finite-volume dependence of these states and 
linking them to the resonances of Nature is a long-term 
program of the lattice QCD community. 

In this letter we report the excited-state energy spec- 
trum of the nucleon in the light quark-mass regime of 
QCD for the first time. Of particular note is the iden- 
tification of a new low-lying state associated with the 
Roper resonance of Nature. 

Several attempts have been made in the past to find 
the elusive low-lying Roper state in the lattice frame- 
work |I3-22l]. The results reported therein are as much 
about the development of lattice techniques as they are 
about the nucleon spectrum. Where the lattice tech- 
niques are regarded as robust, a low-lying Roper state 
was not observed. The difficulties lie in finding effec- 
tive methods to isolate the energy eigenstates of QCD 
and in accessing the light quark mass regime of QCD. 

The 'Variational method' J23, 24] is the state-of-the- 
art approach for determining the excited state hadron 
spectrum. It is based on the creation of a matrix of cor- 
relation functions in which different superpositions of 
excited state contributions are linearly combined to iso- 
late the energy eigenstates. A diversity of excited-state 
superpositions is central to the success of this method. 

Early implementations of this method, using a vari- 
ety of standard spin-flavor interpolating fields of fixed 
source-distribution size, were not successful in isolat- 
ing energy eigenstates. Instead the putative eigenstates 
were superpositions of energy-states B25I1 and the low- 
lying Roper state was hidden by excited-state contami- 
nations. A solution to this problem was established in 
Refs. 11261 12711 where a low-lying Roper state was iso- 
lated. Key to this approach, used herein, is the uti- 
lization of a diverse range of fermion source and sink 
smearings in creating the matrix of correlation func- 
tions. The diversity of smearings leads to a wide variety 
of superpositions of excited-state contributions, provid- 
ing a suitable basis for constructing linear combinations 
which isolate the eigenstates. 

These effective techniques J25 - 28 ] were developed in 
the quenched approximation and we bring them to the 
dynamics of full QCD for the first time here. The low- 
lying even-parity states of the nucleon are explored in 
full QCD using 2+ 1 -flavor dynamical-QCD gauge-field 
configurations HI] . Whereas other recent full QCD anal- 
yses report a first positive parity excited state 
that appears too high to be considered as the Roper res- 
onance, we will illustrate how the low-lying Roper state 
of the present analysis approaches the physical mass 
of Nature, displaying significant chiral curvature at the 
lightest quark mass (corresponding to a pion mass of 



156 MeV, only slightly above the physical value of 140 
MeV). These full QCD results, providing the world's 
first insight into the nucleon mass spectrum in the light- 
quark regime, are significantly different from those of 
quenched QCD and provide interesting insights into the 
dynamics of QCD. 

In constructing our correlation matrix for the nucleon 
spectrum, we consider the two-point correlation func- 
tion matrix with momentum p = 



Gp) = J] Tr sp {r ± <n^(*)*;(0)|fi>}, 



z 



a iq- -m a t 



A" A 



(i) 



(2) 



where Dirac indices are implicit. Here, A" and A" are the 
couplings of the interpolators xi an d Xj at tne sm k an d 
source respectively and a enumerates the energy eigen- 
states with mass m a . Y ± = i(yo ± 1) projects the parity 
of the eigenstates. 

Since the only f dependence comes from the expo- 
nential term, one can seek a linear superposition of in- 
terpolators, Xj u p sucn tnat 



Gijito+At)^ = e- m ° A 'G i j(t )u«, 



(3) 



for sufficiently large to and fo+Af. Multiplying the above 
equation by [Gy(/o)] _1 from the left leads to an eigen- 
value equation 

[(G(t )r 1 G(t + At)] ij u a j =c a uf, (4) 

where c a = e~ m, ' At is the eigenvalue. Similar to Eq. 
one can also solve the left eigenvalue equation to re- 
cover the v a eigenvector 



v? [G(t +At)m )r 1 ] iJ = c a v*. 



(5) 



The vectors u" and v" diagonalize the correlation matrix 
at time fo and fo + Af and provide the projected correlator 

vfGf/f)/^. (6) 

The parity projected, eigenstate projected correlator 

G a ± =V*Gfft)u a j, (7) 

is then analyzed using standard techniques to obtain the 
masses of different states 

The PACS-CS 2 + 1 flavor dynamical-fermion con- 
figurations [8] made available through the ILDG 
are used herein. These configurations use the 
non-perturbatively (9(a)-improved Wilson fermion ac- 
tion and the Iwasaki-gauge action B30I1 . The lat- 
tice volume is 32 3 x 64, with /3 = 1.90 provid- 
ing a lattice spacing a = 0.0907 fm. Five val- 
ues of the (degenerate) up and down quark masses 
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Figure 1: (Color online). Effective mass from smeared-source to point-sink correlators for various levels of smearings at the source for k u j 
0.13770 (left) and k uc i = 0.13700 (right), corresponding to pion masses of 293 and 702 MeV, respectively. 



are considered, with hopping parameter values of 
K ud = 0.13700, 0.13727,0.13754, 0.13770and 0.13781, 
corresponding to pion masses of m„ = 0.702, 0.572, 
0.413, 0.293, 0.156 GeV fl; for the strange quark k s = 
0.13640. We consider 350 configurations for the four 
heavier quarks, and 198 configurations for the lightest 
quark. An ensemble of 750 samples for the lightest 
quark mass is created by using multiple fermion sources 
on each configuration, spaced to sample approximately 
independent regimes of each configuration. Our error 
analysis is performed using a second-order jackknife 
method, where the ^ 2 /dof for projected correlator fits 
is obtained via a covariance matrix analysis. Our fitting 
method is discussed in Refs. $25, 271. 

The complete set of local interpolating fields for the 
spin- 1 nucleon are considered herein. Three different 
spin-flavor combinations of nucleon interpolators are 
considered, 

X x (JC) = e ahc (u Ta (x)Cy 5 d b (x))u c (x) , (8) 
Xi(x) = ^\u r \x)Cd\x))y s u c {x) , (9) 
X4 (x) = e ahc (u Ta (x)Cy 5 y 4 d h (x))u c (x). (10) 

Each interpolator has a unique Dirac structure giving 
rise to different spin-flavor combinations. Moreover, as 
each spinor has upper and lower components, with the 
lower components containing an implicit derivative, dif- 
ferent combinations of zero and two-derivative interpo- 
lators are provided. The local scalar-diquark nucleon 
interpolator, ^1, is well known to have a good overlap 
with the ground state of the nucleon. Also, this xi in- 
terpolator is able to extract a low-lying Roper state in 
quenched QCD [27[]. On the other hand, the \i interpo- 
lator, which vanishes in the non-relativistic limit, cou- 



ples strongly to higher energy states. The interpolator 
Xa is the time component of the local spin-| isospin-^ 
interpolator which also couples to spin-i states. It pro- 
vides a different linear combination of zero- and two- 
derivative terms. 

In constructing our correlation matrices, we first con- 
sider an extensive sample of different levels of gauge- 
invariant Gaussian smearing i3ltl at the fermion source 
and sink, including 4, 9, 16, 25, 35, 50, 70, 100, 125, 
200, 400, 800 and 1600 sweeps. These levels of smear- 
ing correspond to rms radii in lattice units (a 0.09 fm) 
of 1.20, 1.79, 2.37, 2.96, 3.50, 4.19, 4.95, 5.920, 6.63, 
8.55, 12.67, 15.47 and 16.00. 

FigQ] displays effective mass plots, m(t) = 
\n{Gt.(t) / Gt(t +1)} for smeared-source to point-sink 
correlators of X \- The variation in the superposition 
of excited state contributions is revealed in the differ- 
ent approaches of the effective mass to the ground state 
plateau where the results converge. From these plots, it 
is clear that the correlation matrix analysis for excited 
state contributions will be most effective in the regime 
t = 17 - 21 where there is diversity in the curves. 

To explore the low-lying eigenstates of the nucleon 
spectrum, we construct several 4x4 correlation matrices 
of^fi as described in TableQ] These matrices provide ro- 
bust results for the lowest three energy eigenstates ob- 
served, with the highest energy level accommodating 
the fourth eigen-energy and any residual strength from 
higher states not eliminated via Euclidean time evolu- 
tion. 

The masses from the projected correlation functions 
obtained from the correlation matrix analysis are very 
consistent over the variational parameters (fo,Af) as il- 
lustrated in Fig. |2] Careful examination of Fig. [2]re- 
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Table 1: Smearing levels used in constructing 4x4 correlation matrix bases. 
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veals some systematic drift in the second excited state 
mass at small At for to = 17, 18. This emphasises the 
preference for sele cting l arger values of (to, At) as dis- 



selecting lar 

cussed in Refs. However, larger uncertain- 

ties are evident for to = 18, 19 with large At due to the 
suppression of excited states via Euclidean time evolu- 
tion. We select fo=18, At-2 as providing the best bal- 
ance between these systematic and statistical uncertain- 
ties M25L 12711 . These variational parameters also provide 
projected correlation functions having favorable ;f 2 /dof 
in the effective mass fits. 

The consistency of the extracted masses from all the 
4x4 matrices considered in Table [TJ is illustrated in 
Fig. [3] In particular, the ground and Roper states are 
robust. Both lower and higher smearing radii are bene- 
ficial for spanning the space of states at all quark mass. 
However, we avoid bases which include extreme smear- 
ing counts (400 and 800) as these often provide ill— 



> 

02 



o g.s 
a l al e.s 
□ 2 nd e.s 



5 5 5 5 555? I i I 

OOOO DOOO ooo 



17 



18 



19 



20 



Figure 2: (Color online). Masses from the projected correlation func- 
tions as shown in Eq.0 for each set of variational parameters to (major 
axis) and Af (minor axis). This figure corresponds to the lightest quark 
mass, k u j = 0.13781, for which m n = 156 MeV, and the 3rd basis of 
Tabled 



defined correlation matrices. Hence, we select basis 
number 3 as the focus of subsequent analysis. 

In Fig. [4] results for the first 12 eigenstate energies are 
reported for the five quark masses available. The scale 
is set via the Sommer parameter |33|]. The 12 states 
are drawn from three 8x8 correlation-matrix analyses 
for pairs of xu Xi an d X4- The matrices are formed 
with each interpolator having four levels of smearing. 
Whereas the x\ > X2 and X4, Xi analyses reveal the same 
spectrum, four new states are revealed in the^i,^4 anal- 
ysis in place of the four states dominated by xi- 

The colour coding and symbols for the states of Fig. [4] 
illustrate the flow of the states from one quark mass 
to the next as identified by the eigenvector annihilat- 
ing the state. The structure of the eigenvectors isolat- 
ing the states varies only slowly from one quark mass 
to the next making it easy to trace the propagation of 
the states from the heavy to the light quark-mass region. 
We note that this analysis reveals the same number of 
energy states as in the physical resonance spectrum be- 
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Figure 3: (Color online). Masses of the 2Vi energy states for various 
4x4 correlation matrix bases as given in Tablefl] for Kud = 0.13770, 
m n = 293 MeV, over 50 configurations. 
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Figure 4: (Color online) Masses of the positive-parity states of the nucleon. Physical values are obtained from Ref. [32] and plotted at the far left. 
We note the presence of three approximately degenerate states for the Roper at the third quark mass which split at the lighter quark masses. Lattice 
results for the Roper (red triangles) reveal significant chiral curvature towards the physical mass. Non-interacting P-wave and S -wave threshold 
states are also presented. 



low 2.2 GeV 132] 



Among the most significant results of this investiga- 
tion is the manner in which the extracted Roper state 
(filled triangles) approaches the physical value. The sig- 
nificant curvature in the chiral regime indicates the im- 
portant role played by mesonic dressings of the Roper. 
Non-interacting P-wave Nn (E^ + E„) and S -wave Nnn 
(MN+M n +M n ) threshold scattering states are presented 
by the dashed and dotted lines, respectively. For the 
two large quark masses, the results sit close to the P- 
wave Nn scattering threshold whereas the masses for 
the lighter three quark masses sit much higher. While 
there is some evidence for the P- and S -wave scattering 
states for the two heavier quarks, there is no evidence of 
these states at light quark masses. 

A possible explanation for this feature is that the 
attractive mass-dependent and spin-dependent forces 
which are necessary for the formation of a strong 
resonance only have sufficient strength at light quark 
masses. For example, it is typical to encounter spin- 
dependent forces which are inversely proportional to 
the product of the quark masses undergoing gluon 
exchange. As there is no evidence for the thresh- 
old scattering states with back-to-back momenta of 



one lattice unit, p = (2n/L x ,0,0), at the lightest 
three quark masses, we find it unlikely that scatter- 
ing states with the next back-to-back momenta, p = 
(2n/L x , 2n/L y , 0) would suddenly appear in our spec- 
trum. At light quark masses, resonant eigenstates dom- 
inated by single-particle states dominate the spectral 
function whereas at heavy quark masses, only the multi- 
particle states have spectral strength sufficient to be seen 
in the spectrum. 

In addition to this quark-mass effect, there is a well- 
known volume effect. The couplings to the multi- 
particle meson-baryon states are suppressed by 1 / -VV 
relative to states dominated by a single-particle state. 
On our relatively large volume, it is likely that multi- 
particle states will be suppressed and missed in our 
spectrum, particularly at lighter quark masses where 
the quark-mass effect also acts to suppress the spectral 
strength. Further analysis of finite volume effects (34 ] 
on the spectrum is desirable. 

Future calculations should also investigate the use of 
five-quark operators to ensure better overlap with the 
multi-particle states. This type of novel work using the 
stochastic LapH method is in progress 13511 . Indeed, 
complete knowledge of the spectrum is required for a 
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Figure 5: (Color online). A comparison of the low-lying positive- 
parity spectrum of dynamical QCD (full symbols) and quenched QCD 
results (open symbols) from Ref. 12611 . In the quenched case, the lat- 
tice volume was (2 fm) 3 . 

definitive determination of the properties of the Roper 
resonance. 

Fig. [5] provides a comparison of our results in full 
QCD with earlier results in quenched QCD |2ql . where 
the effects of dynamical quark loops are not considered. 
While the ground state of the nucleon in quenched (open 
symbols) and full QCD (full symbols) are in reason- 
able agreement, significant differences are observed for 
the Roper in the light quark mass regime. This is an 
important discovery emphasizing the role of dynamical 
fermion loops in the structure of the Roper. Note, the 
difference between the f-wave nN threshold scattering 
state energies in full and quenched QCD is due to the 
difference in the lattice volume. 

We have also compared [36] our extracted nucleon 
spectrum with the HSC collaboration's 113711 results at 
m n 2 x, 0.27 GeV 2 where multi-particle states are seen in 
our analysis. Taking into account their small-volume of 
(1.97 fm) 3 , which shifts the energy of the f-wave scat- 
tering states for example, both spectra are in good quali- 
tative agreement. Of particular note is the identification 
of the same number of energy states below 3 GeV. 

This investigation is the first to illustrate the man- 
ner in which the Roper resonance of Nature manifests 
itself in today's best numerical simulations of QCD. 
The quark mass dependence of the state revealed herein 
substantiates the essential role of dynamical fermions 
and their associated non-trivial light-mesonic dress- 
ings of baryons, which give rise to significant chiral 
non-analytic curvature in the Roper mass in the chiral 
regime. 
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